1/1 


IID-fll67  211 
UNCLASSIFIED 


PROBLEHS  IN  NONLINEAR  ACOUSTICS:  PARANETRIC  RECEIVING 
ARRAYS  FOCUSED  FINI. .  (U>  TEXAS  UNIV  AT  AUSTIN  DEPT  OF 
HECHANICAL  ENOINEERINQ  N  F  HAHILTON  20  APR  8C 
N0M14-8S-K-G7G8  F/G  28/1 


PROBLEMS  IN  NONLINEAR  ACOUSTICS: 


PARAMETRIC  RECEIVING  ARRAYS. 
FOCUSED  FINITE  AMPLITUDE  SOUND. 

AND  DISPERSIVE  NONLINEAR  INTERACTIONS 

First  Annual  Summary  Report  under  ONR 
Contract  N00014-85-K-0708  (April.  1986) 

Mark  F.  Hamilton 


S 


DTIC 

ELECTE 
MAY  1  3 1986 


DISTRIBUTION  STATEMENT  K 


Approved  tot  public  tel»as«| 
Diftctbutlon  Uolimitod 


%  S'  /A  /6/ 


UNCLASSIFIED 


SECuniTy  classification  of  this  page  CltTiin  0*la  Bnfnd) 


REPORT  DOCUMENTATION  PAGE 

READ  INSTRL'CTIO."(S 

BEFORE  COMPLETING  FORM 

s.  RECIPIENT'S  catalog  NUMBER 

J _ 

4.  title  fAnd  Sub(i»/»j  Problems  in  Nonlinear  Acoustics: 
Parametric  Receiving  Arrays,  Focused  Finite  Am¬ 
plitude  Sound,  and  Dispersive  Nonlinear  Inter¬ 
actions.  First  Annual  Summary  Report  under 
Contract  N00014-85-K-0708 

S.  TYPE  OF  REPORT  A  PERIOD  COVERED 

Annual  Summary  Report 

1  Sept.  1985  -  30  April  1986 

§■  PERFORMING  ORG.  REPORT  NUMBER 

7.  author^*; 

Mark  F.  Hamilton 

B.  contract  or  grant  NUMBERfaJ 

N00014-85-K-0708 

s.  performing  organization  name  ano  address 

Department  of  Mechanical  Engineering 

The  University  of  Texas  at  Austin 

Austin,  TX  78712-1063 

10.  program  ELEMENT,  project.  TASK 
AREA  A  WORK  UNIT  NUMBERS 

61153N  11 

4126943 

11.  CONTROLLING  OFFICE  NAME  ANO  ADDRESS 

Office  of  Naval  Research 

Department  of  the  Navy 

Arlington,  Virginia  22217 

12.  REPORT  DATE 

30  April  1986 

IS.  NUMBER  OF  PACES 

29 

14.  monitoring  agency  name  a  ADORESSri/  dllUnnt  from  Conltolllng  Officm) 

IS.  SECURITY  CLASS,  (ot  thim  tmpori) 

UNCLASSIFIED 

ISa.  DECLASSIFtr  ATION'  OOWNCRAOlNC 

schedule 

1«.  OlSTRieuTlON  STATEMENT  (ot  (hit  R0pOft) 

Approved  for  public  release;  distribution  unlimited 

(7.  DISTRIBUTION  STATEMENT  ^0/  cbafrccf  wtt«r*d  in  Btock  20,  it  BUtmfmtt  from  RmpoH) 

IB.  supplementary  notes 

IB.  KEY  WORDS  (Coniinxtm  on  rovoroo  mid*  it  nmcmmmmry  m\d  id0nUft  by  BfecA  nu0i6oO 

nonlinear  acoustics 
parametric  array 
focused  sound 
dispersive  interactions 

20.  abstract  (Conllnum  on  ravaraa  gtd9  If  ndcmdmmry  m%d  Idmntity  bjr  AlocA  nuaytmO 

Three  projects  are  discussed  in  this  report.  (1)  Use  of  a  parametric  re¬ 
ceiving  array  to  measure  ship  noise.  The  possibility  of  using  a  parametric 
receiving  array  (parray)  for  eliminating  multipath  components  in  measurements 
of  ship  noise  was  investigated  theoretically.  A  model  based  on  the  non- 
collinear  interaction  of  two  diffracting  Gaussian  beams  has  been  developed. 
Results  were  obtained  for  the  specific  case  of  a  parray  whose  pump  is  located 

(continued  on  back) 

DO  ijShm  1473  edition  OF  1  MOVES  If  OBIOLETE  UNCLASSIFIED 


S/N  0102- LF-  014-  6601 


security  CLAttlFlCATlON  OF  THIS  PAOC  fRIiwi  Of  tnindd) 


UNCLASSIFIED 


SECUNITV  classification  of  this  page  (Whit  0««  Enl*r«4 


on  the  sound  source.  Predicted  parametric  measurements  of  beam  patterns 
radiated  by  the  source  are  shown  to  be  in  reasonable  agreement  with  con¬ 
ventional  results  when  the  receiving  hydrophone  is  beyond  about  10 
Rayleigh  distances  from  the  source.  Inversion  formulas  were  derived  that 
determine  the  noise  source  distribution  from  nearfield  measurements  ob¬ 
tained  with  a  parray.  (2)  Focused  finite  amplitude  sound.  An  existing  com¬ 
puter  program  has  been  modified  for  application  to  focused  sound  fields 
radiated  by  real  sources.  The  program  accounts  simultaneously  for  the 
effects  of  nonlinearity,  diffraction,  and  absorption.  Results  computed 
with  the  program  compare  very  well  with  analytic  linear  models  for 
focused  sound  that  is  strongly  affected  by  diffraction.  Analysis  of  low 
amplitude  focused  sound  fields  was  performed  before  using  the  program  to 
investigate  finite  amplitude  effects  at  high  amplitudes.  (3)  Dispersive 
nonlinear  wave  interactions  in  a  rectangular  duct.  Preliminary  experimen¬ 
tal  work  was  performed  in  preparation  for  an  investigation  of  finite  ampli¬ 
tude  noise  that  propagates  in  higher  order  modes  of  a  rectangular  duct. 


feif 


'^1  irv 

'  ‘  frrto 


M  CIO?-  Lf-  ni.1-  6601 


UNCLASSIFIED 

SECuNITy  classification  of  This  PAOEfW»l»n  Dm*  £nl»r»a> 


INTRODUCTION 


Work  discussed  In  this  Annual  Summary  Report  was  supported  by  ONR 
Contract  N00014-85-K-0708.  The  contract  began  1  September  1985  and 
terminates  31  August  1986,  but  this  report  covers  research  performed 
only  up  to  3D  ^^1986.  The  personnel  Involved  with  the  research  effort 
are  as  follows. 


Senior  Personnel 


M.  F.  Hamilton,  principal  investigator 


J.  Naze  TJetta  and  S.  Tjetta,  visiting  scientists, 
University  of  Bergen,  Norway 


Graduate  Students 

C.  Darvennes.  Ph  D.  student  In  Mechanical  Engineering 
T.  S.  Hart,  MS.  student  in  Electrical  Engineering 
S.  J.  Lind,  MS.  student  In  Architectural  Engineering 


It  must  be  pointed  out  that  the  research  performed  by  J.  Naze  Tjetta 
and  s.  Tjetta  on  problems  outlined  in  the  proposal  for  Contract 
N00014-85-K-0708  was  accomplished  from  June  through  August  of  1985 
and  therefore  Just  prior  to  the  start  of  the  present  contract.  The  TJettas' 
work  was  supported  by  the  Independent  Research  and  Development 
program  of  Applied  Research  Laboratories,  The  University  of  Texas  at 
Austin  (ARL:UT).  Their  collaborative  effort  with  Hamilton  during  that 
period  yielded  substantial  progress  on  the  proposed  use  of  a  parametric 
receiving  array  to  measure  ship  noise.  The  results  of  that  work  are 
therefore  discussed  in  the  present  report. 
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Nr.V*' 


Contract  N000M-85-K-0708  provides  for  basic  research  on 
nonlinear  acoustics.  The  three  projects  currently  supported  by  the 
contract  are  as  follows: 


•*  v\* 


t,  Use  of  a  Parametric  Receiving  Array  to  Measure  Ship  Noise  j 
2,  focused  Finite  Amplitude  Sound  '  ” 

3;  Dispersive  Nonlinear  Wave  Interactions  in  a  Rectangular  Duct . 


The  first  project  is  the  primary  research  effort  supported  by  the 
contract,  and  the  investigation  is  predominantly  theoretical.  The  second 
Involves  computational  analysis,  the  third  experimental.  Each  project  is 
discussed  separately  in  the  following  sections. 
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I.  USE  OF  A  PARAMETRIC  RECEiVINC  ARRAY  TO  MEASURE  SHIP 
NOISE 


1.  Personnel 


Research  on  this  project  began  during  the  months  just  prior  to  the 
start  of  the  present  contract  on  1  September  1985.  Substantial 
contributions  were  made  by  J.  Naze  Tjetta  and  S.  Tjetta  from  June  through 
August  of  1985.  The  Tjettas  were  on  leave  from  the  University  of  Bergen, 
Norway,  and  spent  the  summer  months  at  ARL:UT,  where  they  were 
supported  by  internal  funds.  Together  with  Hamilton,  the  Tjettas 
submitted  for  publication  in  the  Journal  of  the  Acoustical  Society  of 
America  a  manuscript  that  discusses  a  mathematical  model  of  the 
parametric  receiving  array,  for  use  in  measuring  ship  noise,  that  is  based 
on  Gaussian  beam  theory.*  It  seems  appropriate  to  report  here  the 
results  contained  in  that  manuscript.  C.  Darvennes  has  wodced  on  this 
project  since  I  September  1985. 


2.  BacKflcguno 

The  task  of  accurately  measuring  the  radiated  noise  from  large  ships 
underway  is  complicated  by  two  factors^  resulting  from  the  phenomenon 
referred  to  as  multipath.  A  hydrophone  used  to  detect  the  acoustic 
radiation  may  receive  not  only  the  noise  propagating  directly  from  the 
ship  but  also  reflected  signals  from  both  the  surface  of  the  water  and  the 
ocean  floor.  One  problem  arises  when  the  ship  is  underway,  because  the 
Doppler  shift  associated  with  each  multipath  component  will  in  general  be 
different.^  The  resulting  differential  Doppler  shift  alters  the  frequency 
distribution  of  acoustic  energy  in  the  received  signal.  To  avoid 
differential  Doppler  shift,  a  highly  directional  hydrophone  array  capable 
of  discriminating  against  the  multipath  components  is  desirable.  The 
second  complication  resulting  from  multipath  concerns  the  determination 
of  transmission  loss  necessary  for  calculating  source  levels.  In  a  recent 
investigation  by  Wright  and  Cybulski,^  sonobouys  were  used  to  measure 


the  low  frequency  noise  radiated  by  targe  merchant  ships  underway. 
Although  the  water  was  approximately  4  km  deep,  an  inablilty  to 
adequately  account  for  bottom  reflections  significantly  reduced  the 
amomt  of  usable  data.  Because  of  its  relative  insensitivity  to 
differential  Doppler  shift  and  its  potentially  high  directivity,  an  end-fire 
array  presents  itself  as  a  viable  means  of  discriminating  against  the 
undesirable  multipath  components. 

Chotiros^  has  suggested  the  use  of  a  parametric  receiving  array 
(parray)  for  measuring  noise  radiated  by  a  ship  underway.  A  parray 
consists  of  only  two  components.  One  is  a  source,  referred  to  as  the 
pump,  capable  of  radiating  a  collimated,  high  intensity  sound  beam.  The 
other  is  a  hydrophone  located  along  the  acoustic  axis  of  the  pump. 
Because  of  its  high  intensity,  the  pump  wave  interacts  nonlinearly  with, 
and  thus  modulates,  any  noise  signal  propagating  through  its  path.  The 
resulting  intermodulation  components  detected  by  the  hydrophone  provide 
information  about  the  noise  signal.  The  differential  elements  that  form  a 
continuous  end-fire  array  are  thus  synthesized  by  the  nonlinear  response 
of  the  fluid  to  the  pump  wave.  When  the  pump  wave  and  noise  signal  are 
each  a  pure  tone  of  different  frequency,  the  interaction  effects  of 
interest  are  the  sum  and  difference  frequency  signals. 

A  unique  advantage  associated  with  using  a  parray  to  monitor  noise 
radiated  by  a  ship  underway  Is  the  ability  to  attach  one  end  of  the  parray, 
specifically  the  pump,  to  the  ship  itself,  with  the  hydrophone  placed  at  a 
fixed  location  some  distance  away  (see  Fig.  I ).  Navigation  data  provided 
by  the  ship  may  be  employed  to  continuously  direct  the  pump  toward  the 
hydrophone.  In  this  way  the  parray  can  maintain  optimal  orientation  for 
eliminating  multipath  components  from  the  noise  signal  radiated  by  the 
passing  ship.  Other  engineering  details  associated  with  this  parray 
configuration  are  discussed  by  Chotiros.^ 

The  remaining,  unresolved  fundamental  problem  involves  determining 
the  acoustic  radiation  pattern  in  the  farfield  of  the  ship  when  the  parray 
itself  may  be  predominantly  within  the  nearfield.  At  frequencies  below 
100  Hz.  the  nearfield  of  large  ships  can  extend  from  hundreds  to  thousands 
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Low  frequency 
beam  pattern 


Figure  1 

Proposed  configuration  for  a  parametric  receiving  array 
used  to  measure  ship  noise 


I 

of  meters.  An  upper  bound  on  the  length  of  the  parray  is  governed  \ 

primarily  by  the  maximum  range  at  which  the  intensity  of  the  pump  wave 

is  capable  of  maintaining  sufficiently  strong  nonlinear  interaction  with 

the  noise  signal.  Maximum  parray  lengths  are  typically  on  the  order  of 

hundreds  of  meters,  although  design  parameters  for  parrays  tens  of  \ 

kilometers  in  length  have  been  discussed  by  McDonough.^  in  any  event. 

interaction  in  the  nearfield  of  both  the  pump  and  noise  source  must  be 

fully  accounted  for.  Specifically,  methods  must  be  developed  with  which 

the  noise  source  distribution  (and  therefore  the  farfield  directivity,  ^ 

which  is  determined  by  the  Fourier  transform  of  the  source  distribution) 

can  be  determined  from  measurements  in  the  nearfield.  None  of  the 

theoretical  models  of  parrays  previously  reported  in  the  literature 

accounts  for  nearfield  diffraction  of  both  the  pump  wave  and  noise  signal. 

Preliminary  experiments  modeling  the  proposed  application  of  a 
parray  for  measuring  ship  noise  have  been  performed  at  ARL:UT.^  The 
experiments  consisted  of  employing  the  method  described  above  to 
measure  low  frequency  radiation  of  a  pure  tone  from  a  large  circular 
source.  The  small,  high  frequency  pump  was  located  in  the  center  of  the 
source  and  was  rotated  to  scan  the  pump  wave  through  the  entire  sound 
field  (as  in  Fig.  2  but  with  (l  the  displacement  of  the  pump  from  the 
center  of  the  source,  equal  to  zero).  Measurements  of  the  sound  field 
were  made  with  the  hydrophone  both  within  and  beyond  the  nearfield  of 
the  source.  Beam  patterns  for  the  radiated  sound  were  thus  obtained  with 
the  parray  and  compared  with  beam  pattemb  obtained  by  conventional 
measurement  of  the  field.  Excellent  agreement  was  obtained  when  the 
hydrophone  was  located  well  in  the  farfield,  but  as  the  hydrophone  was 
moved  into  the  nearfield  the  beam  patterns  measured  with  the  parray 
increasingly  overestimated  the  beamwidths  measured  by  conventional 
means.  Discrimination  of  the  parray  against  multipath  was  also 
demonstrated  in  the  experiments  by  placing  a  reflector  parallel  to  the 
acoustic  axis  of  the  noise  source.  Although  beam  patterns  measured 
conventionally  in  the  presence  of  the  reflector  were  significantly  altered 
relative  to  the  free  field  patterns,  those  measured  with  the  parray 
closely  resembled  beam  patterns  obtained  in  the  absence  of  the  reflector 
The  experiments  clearly  established  the  viability  of  using  a  parray  to 
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separate  the  beam  patterns  due  to  direct  radiation  from  the  contributions 
due  to  multipath  components. 


3.  Results 

In  what  follows  are  discussions  and  results  taken  directly  from  a 
manuscript  submitted  to  the  Journal  of  the  Acoustical  Society  of  America 
by  Hamilton,  Naze  Tjetta  and  Tjetta.  * 

A  model  was  developed  that  describes  the  noncoil  inear  interaction  in 
a  nondissipative  fluid  of  two  sound  beams  generated  by  Gaussian  shaded 
sources.  The  sources  may  possess  linear  phase  shading  and  different  spot 
sizes,  and  their  centers  may  be  displaced  relative  to  each  other  (see  Fig. 
2).  The  analysis  is  based  on  the  fundamental  theory  developed  in  earlier 
papers  by  the  authors.^*®  Numerical  results  were  obtained  for  the 
specific  case  of  a  parray  whose  pump  is  located  on  the  sound  source.  The 
results  show  the  effect  of  the  pump's  position,  size,  and  frequency  on 
parametric  measurements  when  diffraction  of  both  interacting  waves  is 
important.  Asymptotic  results  were  derived  for  the  case  where  the 
receiving  hydrophone  is  located  in  the  farfield  of  the  source. 

The  geometry  is  as  shown  in  Fig.  2.  Gaussian  shading  of  the  form 
exp{-(|x|/ej)^}  is  assumed  for  the  amplitude  of  both  the  pump  (i-1)  and 

source  (i»2),  where  x  is  a  vector  in  the  plane  of  the  pump.  A  phase  term 
is  introduced  to  steer  the  source  relative  to  the  pump.  Solutions  for  the 
sum  and  difference  frequency  waves  received  at  the  hydrophone  are 
obtained  from  a  parabolic  approximation  of  the  nonlinear  wave 
equation.^'®  The  hydrophone  is  located  on  the  acoustic  axis  of  the  pump. 
Thermoviscous  losses  in  the  fluid  are  ignored,  with  the  result  that  the 
solutions  reduce  to  functions  of  the  exponential  integral  that  are  easily 
evaluated.  In  the  figures  that  follow  It  is  assumed  that  the  radius  of  the 
source  is  20  times  that  of  the  pump,  i.e.  62/6  j=20.  The  size  of  the  pump 

measured  in  wavelengths  of  its  radiated  sound  is  fixed  by  setting  k^e^^lS, 

where  k|  is  the  wavenumber  of  the  pump  wave.  The  frequency  of  the 


source  Is  then  determined  by  the  ratio  Parametrically  measured 

beam  patterns  are  calculated  at  ranges  given  by  (5>^il22,  where 
Z2*k2®2^^2  is  the  collimation  length  of  the  beam  formed  by  the  source 

waves.  The  transition  region  between  the  near  and  far  fields  of  the 
source  is  therefore  located  around  02*  1  ■ 

In  Fig.  3  the  solid  curves  are  numerical  results  for  beam  patterns  of 
the  difference  frequency  wave,  referred  to  here  as  parametric  beam 
patterns,  measured  by  the  hydrophone  when  the  pump  is  located  in  the 
center  of  the  source.  The  parametric  beam  patterns  are  calculated  for  O2 

equal  to  1.  10,  and  100,  and  for  source  frequencies  determined  by 

equal  to  3,  10,  and  30.  The  dashed  lines  in  Fig.  3  are  the  linear  farfleld 
directivity  patterns  for  the  source  waves,  which  are  given  by 
exp{-(k2e2Sina)2/4).  It  is  seen  that  for  each  source  frequency,  the 

parametric  pattern  approaches  the  desired  linear  farfield  pattern  as  the 
hydrophone  is  moved  farther  from  the  source.  Beyond  02*10,  the 

parametric  patterns  are  within  1  dB  of  the  linear  farfield  patterns 
between  the  -3  dB  points  of  the  latter.  In  all  cases  the  parametric 
patterns  are  broader  than  the  linear  farfield  patterns,  in  agreement  with 
experiment.® 

The  effect  of  displacing  the  pump  relative  to  the  center  of  the  source 
is  shown  in  Fig.  4.  The  pump  is  located  at  the  edge  of  the  source  (i.e.,  at 
1x1-62),  translation  being  in  the  same  plane  as  that  in  which  the 

source  is  rotated.  It  is  seen  that  displacement  of  the  pump  significantly 
affects  the  parametric  beam  patterns  only  when  the  hydrophone  is  on  the 
order  of  one  Rayleigh  distance  from  the  source.  Other  results,  not  shown 
here,  demonstrate  that  focusing  of  the  pump  wave  on  the  receiver  and 
variation  in  the  size  of  the  pump  have  no  appreciable  effect  on  the 
measured  beam  patterns. 

Finally,  an  inversion  formula  was  derived  that  permits  evaluation  of 
the  source  function  if  the  difference  frequency  pressure  is  known  at  a 
sufficient  number  of  points  at  any  given  range.  The  inversion  formula  is 
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essentially  the  spatial  Fourier  transform  of  the  parametrically  measured 
beam  pattern,  with  a  normalization  factor  that  depends  on  range, 
frequency,  and  the  shape  of  the  pump.  Once  the  source  function  Is 
determined,  the  linear  farfield  directivity  of  the  source  Is  given  simply 
by  a  second  Fourier  transform.  In  theory,  therefore,  parametric 
measurements  made  In  the  nearfield  of  the  source  can  be  used  to 
determine  the  desired  farfield  directivity  pattern.  The  distance  between 
the  pump  and  hydrophone  need  only  be  so  long  as  to  eliminate  unwanted 
multipath  components.  However,  difference  frequency  measurements 
must  be  made  over  a  surface,  not  Just  a  line,  and  both  amplitude  and  phase 
must  be  measured. 

On  1  September  1985,  C.  Darvennes  began  research  on  the  use  of  a 
parray  near  reflecting  surfaces.  One  envisioned  use  of  the  parray  Is  for 
the  measurement  of  noise  radiated  by  surface  ships.  It  Is  therefore 
Important  to  understand  the  extent  to  which  sound  reflected  from  the 
surface  of  the  water  affects  the  response  of  the  parray.  Her  analysis 
begins  with  research  by  Novikov  et  al.^  on  the  use  of  a  parametric 
transmitting  array  located  adjacent  to  a  reflecting  surface.  She  Is  now  in 
the  process  of  verifying  their  analysis  and  attempting  to  modify  the 
results  for  application  to  parametric  receiving  arrays. 
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II.  FOCUSED  FINITE  AMPLITUDE  SOUND 


1.  Personnel 

Work  on  this  project  was  performed  by  T.  S.  Hart  since  1  September 
1985. 


2.  BacKground 

In  a  recent  investigation  of  acoustic  microscopy,  where  focused 
sound  is  used  for  imaging,  questions  were  raised  concerning  the  improved 
resolution  achieved  at  finite  amplitudes.^^'^  ^  Rugar^^  observed  that  at 
sufficiently  high  sound  power  levels,  resolution  beyond  the  linear 
diffraction  limit  can  be  attained.  His  conclusion  is  that  the  nonlinear 
generation  of  higher  harmonics  provides  signals  with  shorter 
wavelengths  than  are  originally  transmitted.  Despite  the  fact  that  it  is 
the  fundamental  transmitted  frequency  that  is  monitored  and  not  the 
higher  harmonics,  improved  resolution  occurs  nevertheless.  Rugar 
suggests  that  after  sound  passes  through  a  focal  region,  the  resulting 
phase  shifts  cause  energy  in  higher  harmonics  to  be  returned  to  the 
fundamental.  As  a  result,  the  information  content  in  the  higher  harmonic 
waves  is  ultimately  imparted  to  the  wave  at  the  fundamental  frequency. 
Rugar's  theoretical  investigation  of  this  phenomenon  is  based  on  a 
quasi  linear  analysis  of  the  second  harmonic  component  formed  by  a 
focused  Gaussian  beam,  conclusions  drawn  from  the  analysis  are 
therefore  necessarily  incomplete.  First,  a  quasilinear  analysis  can  lead 
only  to  indirect  information  about  the  fundamental  component.  Second, 
the  assumption  of  a  Gaussian  source  excitation  obscures  the  intricate 
effects  of  diffraction  that  critically  affect  the  nonlinear  distortion  of 
real  sound  beams. 

A  review  of  published  work  on  focused  finite  amplitude  sound  reveals 
that  the  effects  of  nonlinearity  and  diffraction  have  not  been  considered 
simultaneously  when  both  effects  are  strong.  Within  the  focal  region. 
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however,  both  effects  can  be  quite  significant.  An  early  analysis  of 
spherical  finite  amplitude  waves  by  Naugornykh.  Soluyan.  and  Khokhlov 
was  compared  with  experiment  by  Smith  and  Beyer.  The  spherical  cap 
radiator  used  in  the  experiments  gives  rise  to  strong  diffraction  effects, 
namely,  a  highly  oscillatory  axial  field.  Because  the  spherical  wave 
theory  does  not  account  for  diffraction,  it  is  shown  to  provide  only  crude 
agreement  with  experiment  in  the  prefocal  region,  and  it  is  of  no  use 
within  the  focal  region.  Ostrovskii  and  Sutin*^*'^use  a  spherical  wave 
model  to  describe  propagation  up  to  Just  before  the  focal  region.  Within 
the  focal  region  they  neglect  nonlinearity  and  consider  only  diffraction 
effects.  Consequently,  the  analysis  does  not  properly  account  for  the 
nonlinear  losses  in  the  focal  region  discussed  previously  by  Naugol'nykh 
and  Romanenko. 

Quasi  I  inear  models  have  been  developed  for  systems  where 
diffraction  effects  are  strong  but  finite  amplitude  effects  are  weak.  The 
linear  (fundamental)  and  quasilinear  (second  harmonic)  fields  have  been 
investigated  both  theoretically  and  experimentally  by  Lucas  and 
f^lr  17,18  theory  is  based  on  a  parabolic  approximation  of  the 
nonlinear  wave  equation,  and  agreement  with  experiment  is  very  good. 
Similar  quasilinear  theories  have  been  developed  for  focused  sum  and 
difference  frequency  sound.  ^^”2' 

Only  Bakhvalov,  Zhileikin,  and  Zabolotskaya^^  account  simultaneously 
for  diffraction  and  nonlinearity  in  focusing  systems  where  amplitudes  are 
sufficiently  high  that  the  quasilinear  models  are  inapplicable.  Their 
analysis  is  based  on  numerical  solutions  of  the  parabolic  nonlinear  wave 
equation,  and  it  is  apparently  valid  for  arbitrary  source  excitations. 
Curiously,  however,  the  authors  investigate  only  systems  where 
diffraction  is  of  moderate  importance.  Specifically,  cases  involving 
highly  oscillatoru  sound  fields  that  arise  in  practical 
applications  ^  *  7,  IB  gpg  considered. 


3.  BgSUltS 


A  numerical  Investigation  of  focused  finite  amplitude  sound  based  on 
Aanonsen’s  finite  difference  solutlon^^  of  the  nonlinear  parabolic  wave 
equation  Is  currently  In  progress.  Aanonsen's  algorithm  accounts 
simultaneously  for  nonlinearity,  diffraction,  and  absorption,  but  practical 
considerations  of  computation  time  limit  its  use  to  the  nearfleld  of 
directive  sound  sources.^^  However,  an  appropriate  transformation  was 
subsequently  applled^^  that  facilitates  computations  In  the  farfleld.  and 
excellent  agreement  was  found  between  theory  and  experiment.  The 
transformation  Introduces  a  coordinate  system  that  follows  the  eventual 
divergence  of  sound  beams  In  the  farfleld.  In  the  present  investigation,  a 
similar  transformation  Is  Introduced  that  is  more  suitable  for  describing 
the  convergent  geometry  of  a  focused  system. 

The  geometry  of  the  system  Is  shown  in  Fig.  5.  where  d  Is  the 
location  of  the  focal  plane  along  the  z  axis.  The  system  is  uniquely 
defined  by  three  dimensionless  parameters,  the  linear  focusing  gain  G. 
the  source  strength  B,  and  the  thermoviscous  absorption  A.  in  what 
follows,  thermoviscous  absorption  Is  Ignored.  The  linear  focusing  gain  G 
Is  given  by* ^G-Zq/O.  where  Zq  Is  the  Rayleigh  distance  of  the  source. 

We  define  the  source  strength  B-d/Zg  as  the  ratio  of  the  focal  length  to 
the  plane  wave  shock  formation  distance  Zg-l/pek.  where  p  Is  the 

coefficient  of  nonlinearity,  e  the  acoustic  Mach  number  at  the  source, 
and  k  the  wavenumber.  Ail  pressure  amplitudes  are  normalized  to  the 
amplitude  of  the  fundamental  on  the  source.  For  a  linear  system  with  no 
thermoviscous  losses,  the  behavior  of  the  radiated  sound  field  is 
determined  solely  by  the  gain  G. 

We  first  consider  a  source  whose  amplitude  Is  uniform  within  a 
circle  of  radius  a  and  zero  elsewhere.  In  Fig.  6  the  solid  lines  are  linear 
solutions  obtained  numerically  for  par^eters  similar  to  those  of  the 
experiments  performed  by  Lucas  and  Muir,  *  *  ®l.e.,  G-40.3  and  B-0. 1 1 .  The 
broken  lines  are  the  corresponding  analytical  solutions  of  the  parabolic 
wave  equation.  The  top  curve  is  the  amplitude  of  the  fundamental 
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component  along  the  z  axis,  and  the  middle  curve  Is  the  corresponding 
phase.  Agreement  between  the  numerical  and  analytical  solutions  is  seen 
to  be  very  good  In  the  focal  region,  where  the  amplitudes  differ  by  less 
than  0.2  dB.  In  the  bottom  curve  Is  shown  the  beam  pattern  in  the  focal 
plane,  where  again  the  numerical  and  analytical  solutions  are  very  close 
(the  label  on  the  abscissa  should  read  26(r/a),  where  r  Is  the  radial 
coordinate  perpendicular  to  the  z  axis).  Of  particular  significance  is  the 
fact  that  the  analytical  solutions  are  shown  by  Lucas  and  Mulr^^to  be  In 
very  good  agreement  with  experiment. 

In  Fig.  7  are  quaslllnear  solutions  obtained  numerically  for  the 
second  harmonic  component  that  Is  generated  by  the  sound  field  modeled 
In  Fig.  6.  Whereas  the  maximum  normalized  amplitude  of  the  fundamental 
component  depends  only  on  (and  Is  approximately  equal  to)  6,  that  of  the 
second  harmonic  also  depends  on  the  amplitude  of  the  fundamental.  In  this 
case  the  second  harmonic  signal  attains  an  amplitude  of  about  13  times 
that  of  the  fundamental  component  at  the  source.  In  addition,  the  beam 
pattern  In  the  focal  plane  Is  significantly  narrower  than  that  of  the 
fundamental. 

Finally.  In  Figs.  8  and  9  are  shown  results  for  the  fundamental  and 
second  harmonic  component  generated  by  a  focused  source  whose 
amplitude  varies  as  [l-(r/a)^H  when  ria,  and  which  Is  zero  when  r>a. 
This  source  condition  was  used  by  Baklwalov  et  al.^^  to  Investigate 
diffraction  In  their  numerical  analysis  of  focused  finite  amplitude  sound. 
The  values  of  6  and  B  are  the  same  as  those  used  in  Figs.  6  and  7. 
Comparison  with  Figs.  6  and  7  reveals  that  the  highly  oscillatory  field 
structure  that  characterizes  real  uniform  circular  sources  Is  all  but 
absent  In  the  present  model.  Although  the  effects  of  diffraction  In  the 
focal  region  Increase  with  the  linear  focusing  gain,  the  parameter  G  used 
here  Is  more  than  5  times  the  maximum  value  considered  by  Bakhvalov  et 
al.  As  a  result,  continuation  of  the  present  Investigation  would  yield  the 
first  theoretical  results  for  focused  sound  beams  where  the  effects  of 
both  nonlinearity  and  diffraction  are  strong. 
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Figure  8 

Linear  sound  field  generated  by  a  focused  source  whose  amplitude 
varies  as  (l-ir/a)^]^  when  ria  and  which  is  zero  when  r>a. 
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Figure  9 

Quasilinear  second  harmonic  sound  field  generated 
by  the  linear  sound  field  shown  in  Fig.  8. 
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III.  DISPERSIVE  NONLINEAR  WAVE  INTERACTIONS  IN 
RECTAN6ULAR  DUCT 


t.  Personnel 

This  research  is  to  be  an  extension  of  work  currently  supported  by 
ONR  Contract  N00014-84-K-0574.  As  discussed  by  D.  T.  Blackstock^® 
principal  investigator  under  that  contract,  the  project  currently  involves 
an  investigation  of  finite  amplitude  pure  tones  that  interact  in  higher 
order  modes  of  a  rectangular  duct.  Since  I  January  1986,  S.  J.  Lind  has 
assisted  with  the  completion  of  that  work.  Beginning  1  June  1986,  Lind 
will  extend  the  research  to  include  bands  of  random  noise  that  interact  in 
higher  order  modes. 


^  2.  Ba£isgr.Qund 

Guided  acoustic  wave  propagation  is  of  interest  in  a  wide  variety  of 
physical  systems  that  involve  the  transmission  of  sound  in  bounded  or 
layered  media.  Of  the  many  examples  that  may  be  cited  include  the 
m  propagation  of  sound  in  pipes,  ducts,  shallow  water  channels,  and  surface 

acoustic  wave  devices.  Inherent  to  all  waveguide  systems  is  the  property 
of  dispersion,  the  phenomenon  whereby  the  phase  speed  of  a  wave  depends 
>  on  frequency.  At  the  cutoff  frequency,  which  is  the  lowest  frequency  at 

which  sound  can  propagate  in  any  given  mode,  the  phase  speed  of  sound 
r  down  the  waveguide  is  infinite.  As  the  frequency  is  increased,  the  phase 

speed  decreases  asymptotically  toward  the  speed  of  sound  in  unbounded 
media.  Waveguides  therefore  constitute  systems  that  are  very  highly 
dispersive. 

/  Whereas  the  propagation  of  low  amplitude  sound  in  higher  order 

modes  of  waveguides  is  well  understood,  that  of  high  amplitude  sound  is 
not.  Specifically,  dispersion  induced  by  the  boundaries  of  a  waveguide 
C  critically  affects  the  distortion  of  finite  amplitude  sound.  Sound  waves 

in  the  lowest  order  mode,  or  plane  wave  mode,  are  practically 


nondlspersive.  and  their  behavior  at  finite  amplitudes  Is  well  understood. 
Although  a  number  of  papers^^"^^  consider  guided  waves  of  finite 
amplitude  In  higher  order  modes,  it  is  significant  that  none  of  the 
analyses  accounts  consistently  for  the  combined  effects  of  nonlinearity, 
dispersion,  and  absorption.  Even  more  remarkable  Is  the  complete 
absence  of  any  published  experimental  work,  particularly  when  one 
considers  the  various  approximations  made  In  the  theoretical  work. 

Substantial  progress,  both  theoretical  and  experimental,  has  been 
made  at  ARLrUT  on  the  problem  of  finite  amplitude  sound  propagation  In 
higher  order  modes  of  ducts.  An  eight  meter  rectangular  duct  was  used 
that  allows  the  simultaneous  generation  of  finite  amplitude  sound  In  the 
plane  wave  mode  and  first  higher  order  mode  (see  Fig.  10).  Two  cases  of 
nonlinear  wave  propagation  have  been  Investigated.^^ 

The  first  case  Involves  the  weak  nonlinear  Interaction  of  a  pure  tone 
In  the  plane  wave  mode  (0,0)  with  a  pure  tone  of  another  frequency  In  the 
first  higher  order  mode  (1,0).  A  quasi  linear  theory^^  was  developed  that 
accurately  predicts  the  amplitudes  of  the  waves  generated  nonllnearly  at 
the  sum  and  difference  frequencies.  Of  particular  significance  Is  that  the 
experiment  demonstrated,  for  the  first  time,  dispersion  Induced  spatial 
oscillations  in  a  nonllnearly  generated  wave  field.  The  oscillations  had 
been  predicted  by  researchers  In  nonlinear  optics,  but  the  shorter  optical 
wavelengths  Involved  prevent  similar  measurements  from  being  made. 

The  second  case  Involves  the  strong  nonlinear  distortion  of  a  single 
pure  tone  in  the  (1,0)  mode.  Three  frequency  ranges  may  be  Identified. 
One.  very  near  cutoff.  Is  characterized  by  standing  waves  formed  across 
the  duct.  Very  little  energy  Is  propagated  down  the  duct.  A  second,  at 
frequencies  far  above  cutoff.  Is  characterized  by  moderate  dispersion. 
All  components  of  the  nonllnearly  distorted  wave  propagate  down  the  duct 
with  approximately  equal  phase  speeds.  Although  a  set  of  coupled 
differential  equations  has  been  derived  for  this  frequency  range.  It  Is  not 
amenable  to  numerical  solution.  The  third  Is  In  a  mid-frequency  range 
where  dispersion  is  sufficiently  strong  that  only  waves  propagating  with 
Identical  phase  speeds  are  strongly  coupled.  In  this  range  considerable 
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Waveguide  apparatus. 


simplification  of  the  set  of  coupled  equations  is  possible,  and  agreement 
between  theory  and  experiment  is  good  all  the  way  through  the  fifth 
nonllnearly  generated  harmonlc.^^ 

Two  observations  may  be  made  regarding  the  research  underway  at 
ARL:UT.  First,  results  from  similar  experiments  have  never  before  been 
published.  Second,  the  previously  published  theories  provide  inadequate 
models  of  the  experiments.  The  analyses  developed  at  ARLiUT,  which 
account  for  the  combined  effects  of  nonlinearity,  dispersion,  and 
absorption,  are  in  a  number  of  cases  In  excellent  agreement  with 
experiment.^® 


3.  BflSUUS 

Since  1  January  1986,  Lind  has  assisted  J.  A.  TenCate,  who  is 
supported  by  ONR  Contract  NOOOM'S^-K'OSy^  In  the  completion  of  work 
on  pure  tones  that  interact  In  higher  order  modes  of  a  rectangular  duct.^® 
Through  collaboration  with  TenCate,  Lind  has  acquired  experience  with 
the  same  experimental  apparatus  that  he  will  use  beginning  I  June  1986  to 
extend  the  work  on  pure  tones  to  bands  of  random  noise. 
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